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Abstract 

We discuss the transition strength between the disordered ground state and the basic low-lying 
triplet excitation for the Heisenberg chain with alternating exchange based on series expansions in 
the interdimer interaction up to 10 th order. We demonstrate that the wave vector dependence of 
the simple dimer approximation is strongly modified in higher orders. 
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Introduction 



Low-dimensional quantum antiferromagnets have received much interest in recent years 
since they serve as model substances allowing to investigate in detail the effects of quantum 
fluctuations and to test theoretical models. 

One important class of materials in this context consists of an assembly of dimers (two 
strongly coupled spins 1/2) which interact sufficiently weak to guarantee that the dimer gap 
does not close. These materials are characterized by a disordered singlet ground state and 
a finite spin gap to triplet excited states. A number of investigations has been performed 
on the simplest model systems of this type, the Heisenberg chain with alternating exchange: 
The dispersion of the lowest singlet-triplet transition as well as the two particle bound states 
have been calculated by series expansions in the interdimer interaction up to 10 th orderBSi'S 
and the corresponding structure factor S(q) up to third order!. Here we supplement the 
analysis of Ref. [1] by presenting cluster expansion results for the dynamic structure factor 
up to 10 th order. 

Since the lowest triplet mode is an isolated excitation, the weight of the corresponding 
S— function, i.e. the contribution of the singlet-triplet transition to the structure factor 5(q) 
is the quantity to be discussed. We will use the following notations: The dynamical structure 
factor for spins localized on a Bravais lattice is defined as 

/oo 
dt e-^(5 Q (q)^(-q)} (1) 
-oo 

where 

S a (q) = ^V 4qR S Q (R) (2) 

R 

is the Fourier transformation of the spin operators at lattice sites R. The superscripts a, (3 
denote the spin components and the brakets (•••) thermal expectation values (which for 
T = reduce to groundstate expectation values (0| • • • |0)). 

If we consider transitions from the groundstate to some isolated eigenstate |n) of Ji with 
energy u n (q), we obtain 5-peaked contributions to the dynamical structure factor 

S^(q,w) = ^(0| l S Q (q)|n)(n|5' 3 (-q)|0)^-^(q)) (3) 

n 

= ^/^(q)%-<q))- (4) 



Inelastic neutron scattering probes directly the transition from the (singlet) groundstate to 
the lowest (triplet) excitation \t) of the quantum disordered compound and we will reduce our 
discussion to this contribution to the dynamical structure factor. Owing to the rotational 
symmetry of the model it is sufficient to calculate /f 2 (q) only and we use the shorthand 
J(q) :=/f(q,l). 

In the following we will give theoretical results for the ID array of interacting dimers. 
This model was treated before!, it is, however, instructive (i) to demonstrate for the ID case, 
that the simple lowest order results are modified by further terms which emerge in higher 
orders and (ii) to present the quantitative changes for the transition strength comparing 
first order to 10 th order results. The hamiltonian of this model is of the following form 

N/2 

7~C = J S2jS2i+l + AS2t+lS2i+2, J > 0. (5) 

i=0 

Here, N denotes the number of spins and periodic boundary conditions are used. In the case 
of A = the groundstate of the system consists of rung singlets on the bonds (2i, 2i+l). Local 
rung singlets can be excited to rung triplets which become gapped dispersive excitations 
when switching on A, < A < 1. In the limit A = 1 we arrive at the well known Heisenberg 
antiferromagnet (HAFM) where the lowest, gapless excitations are pairs of S = 1/2 spinons. 

The one magnon excitation energies u)(q) have be obtained by perturbation expansion 
in A up to 9 th order in Ref. |2| and to 10 th order in Ref. [I] using the cluster expansion approach. 

The dynamical structure factor 

For the chain geometry shown in Fig. p] we get for the intensity up to first order in A the 
following expression 

J(q) = sin 2 ^(l + ±Acos(ga)) +0(\ 2 ). (6) 

Here, q is the projection of the wavevector q in the chain direction. We note that the 
term oc sin 2 (qd/2) is typical for systems consisting of isolated dimers and is known as the 
dimer structure factor. The first order correction in A adds an additional modulation to the 
intensity, which depends on the ratio a = ||d||/||a||. 

Using the cluster expansion method (see Appendix ^) we have systematically calculated 
the series in A for the intensity up to 10 th order. This requires linked clusters consisting of 
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maximum 10 bonds. The resulting series can be split into three different terms: 



/(q) = B c (q, A) + B s (q, A) + A(q, A). (7) 
To illustrate the result we give the series up to fourth orderS 



B c (q, A) = sin 2 ^ ^ \ij cos(jqa) 

4 

B s (q, A) = sin qd Uj sm(jqa) 



(8) 



where 



/i = 1 - —A 2 - -A 3 + —A 4 , m = -\- -A 2 - —A 3 + —A 4 , (9) 

F0 16 32 1536 '^28 192 2304 ' v ; 

^ 2 = —A 2 + —A 3 + -^-A 4 , At3 = lA 3 + i^A 4 , /i4 = ^A 4 , (10) 

P2 16 48 1024 ' ^ 3 64 2304 ' p 1024 ' v ; 

1 7 131 1 95 

^0 = 0, ^ = -A 2 + — A 3 - —A 4 , z, 2 = -A 3 + — A 4 , (11) 

' 1 8 192 4608 96 4608 ' V ; 

93 

u 3 = — A 4 , i/ 4 = 0, (12) 

s 23Q4 4 v ; 



and 



A(q,A) = J-A 4 (cos(qd) -cos(2ga)). (13) 

The terms -B c (q, A) are consistent with previous publications!, whereas B s (q, A) and A(q, A) 
are additional corrections. They originate from a complete expansion of both the ground 
state |0) and the first excited state \t). If one assumes that (0\S 2 i\t) = —(0\S2i+i\t) for the 
matrix elements on even and odd sites, one ends up with -B c (q, A) only. However this is 
only correct to first order and in general we have (0\S2i\t) ^ —(0\S2i+i\t). This is related to 
virtual states which occur during the perturbation expansion and exhibit odd parity under 
the exchange of the triplets. These corrections begin to appear in second order. 

In Figs. H] and [3] we show some typical plots of the intensity /(q) for two different ratios 
a = 10/15, 10/20 and two different coupling strength A = 0.4, 0.8. The difference between 
the zeroth and the first order emerges very clearly. The higher order terms emphasize the 
modulation originating from the two length scales ||a|| and ||d||. 
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Sum rule 



Independently of whether the magnetic system is ordered or not, the total integrated scat- 
tering intensity has a well defined magnitude, determined by the local spin length through 
the following sum rule: 

1=J2 u =S(S + 1). (14) 

J d( i 

In the case of the one dimensional alternating chain the integral reduces to the constant 
part Ho, since only non oscillating terms will survive the integration: 

X=-(l--X 2 --X 3 + —X A + ..)<- S = -, A«l. (15) 
4 V 16 32 1536 J ~ 4 2' K J 

For the noninteracting case (A = 0) the sum rule is exhausted by the one triplet excitation 
since it is an exact eigenstate: From (H) we see that this excitation gives the only 
non- vanishing matrix element. Switching on the coupling between the dimers, more and 
more intensity goes in two or more magnon scattering processes. 

Conclusion 

Using the S = \ Heisenberg chain with alternating exchange interactions we have demon- 
strated by high order series expansions that the simple dimer structure factor is considerably 
modified in higher orders in the interdimer exchange. In a forthcoming publication we will 
extend our approach to cover interacting dimers in a threedimensional lattice in order to 
allow a discussion of the transition strength for the interacting dimer materials KCUCI3 
and TICUCI3. These materials have been investigated in detail in the last years in inelas- 
tic neutron scattering (INS) experimentsi'i, thus directly exploring the basic singlet-triplet 
transition in all of reciprocal space. 

APPENDIX A: CLUSTER EXPANSION 

In this appendix we briefly summarize the method of cluster expansion for the dynamical 
structure factor in the case of the alternating chain. Some detailed considerations can be 
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found elsewhere!. As shown in FIG. |l]the cristallographic unit cell contains two spins. Thus 
the Fourier transform of the spin operator splits into two parts and reads as: 

S a (q) = e~ iqna (e- iqd ^(n) + e iqd ^(n)) . (Al) 

n 

As before ||a|| = a is the distance between neighbouring spins and 2d denotes the spacing 
between the two spins on a dimer. In our notation q is the projection of the wave vector q 
on the chain direction. 

Using translational invariance with respect to the center of the dimer we obtain for the 
singlet-triplet transition amplitude: 

/(q) = e~ igan UnH + A Z 2 Z 2 (n) + e 2iqd A^(n) + e~ 2iqd A a 2 { '(n)l (A2) 

n 

where 

A$(n) = (0\S?(0)\1)(1\S? (n)|0>, i,j = 1,2. (A3) 

Here, the sum is taken over all integer numbers. However, it is more convenient to calculate 
the functions Af z (n) for positive numbers n. This is feasible by the benefit of inversion 
symmetry concerning a dimer center of mass which is manifest in the relations: 

AH(-n) = A z 2 z 2 (n) and Ag(-n) = A%(n). (A4) 
Inserting ( |A4| ) into ( |A2| ) one arrives at the following result: 
7(q) = 2 (^u(n) + ^SW) cos(gna) 

n>0 

+ 2 cos(qna - 2qd) + Aff(n) cos(qna + 2qd)) (A5) 



n>0 



+ Aff(O) + A Z 2 Z 2 (0) + e 2iqd A^(0) + e~ 2iqd ^(0). 

We state that the functions A z ?(n) have to be calculated for positive n only. 

At first glance the functions A zz (n) are groundstate expectation values which can be 
computed by well known cluster expansion methodei. The projection operator V = |1)(1| 
has to be evaluated from the one magnon states \ip®) where i labels the lattice site. By 
means of degenerate cluster expansion these states are generated order by order!. Then we 
find for projection operator: 

P=|1><1| = E^" 1 )^ W >^ W I- 



g is the overlapping matrix of the \ip^>): 



(A7) 



To invert g we use the fact that g is the unit matrix for A — > 0: 



9 = 1 + 9- 



(A8) 



Owing to the matrix norm \\g\\ < 1 we apply a geometric series to invert g: 



oo 



(A9) 



Now we have everything at hand to calculate the singlet-triplet-intensity of the dynamical 
structure factor: Apply degenerate perturbation theory to obtain the states \if)^') and V. 
Then calculate g and invert this matrix by using (|A9|) . Finally, apply non degenerate 
perturbation theory to compute the functions Af?(n). 
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FIG. 3: Structure factor with A = 0.4, 0.8 and ratio d/a = 10/15 
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